home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / TUTEST.PAS < prev    next >
Pascal/Delphi Source File  |  1991-05-01  |  791b  |  24 lines

  1. PROCEDURE tutest(data1: glnarray; n1: integer; data2: glmarray;
  2.         n2: integer; VAR t,prob: real);
  3. (* Programs using routine TUTEST must define the types
  4. TYPE
  5.    glnarray = ARRAY [1..n1] OF real;
  6.    glmarray = ARRAY [1..n2] OF real;
  7. in the main routine, with n2 less than or equal to n1.   *)
  8. VAR
  9.    var2,var1,df,ave2,ave1: real;
  10.    i: integer;
  11. BEGIN
  12.    IF (n2 > n1) THEN BEGIN
  13.       writeln('pause in routine TUTEST');
  14.       writeln('first array must be the larger'); readln
  15.    END;
  16.    avevar(data1,n1,ave1,var1);
  17.    FOR i := 1 TO n2 DO data1[i] := data2[i];
  18.    avevar(data1,n2,ave2,var2);
  19.    t := (ave1-ave2)/sqrt(var1/n1+var2/n2);
  20.    df := sqr(var1/n1+var2/n2)/
  21.       (sqr(var1/n1)/(n1-1)+sqr(var2/n2)/(n2-1));
  22.    prob := betai(0.5*df,0.5,df/(df+sqr(t)))
  23. END;
  24.